Efficiency and Its Bounds for Thermal Engines at Maximum Power using a Newton's Law of 

Cooling 
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We study a thermal engine model for which Newton's cooling law is obeyed during heat transfer 
processes. The thermal efficiency and its bounds at maximum output power are derived and dis- 
cussed. This model, though quite simple, can be applied not only to Carnot engines but also to four 
other types of engines. For the long thermal contact time limit, new bounds, tighter than what were 
known before, are obtained. In this case,this model can simulate Otto, Joule Brayton, Diesel, and 
Atkinson engines. While in the short contact time limit, which corresponds to the Carnot cycle, the 
same efficiency bounds as Esposito et al's are derived. In both cases, the thermal efficiency decreases 
as the ratio between the heat capacities of the working medium during heating and cooling stages 
increases. This might provide instructions for designing real engines. 

PACS numbers: 05.70.Ln,05.20.-y 



INTRODUCTION 



It is well known that real thermal engines can not 
achieve a perfect Carnot cycle. In a perfect Carnot 
cycle,the two reversible isothermal stages must be in- 
finitely long and hence the Carnot engine has zero 
power output. Although the Carnot thermal machine 
is impractical, it gives an upper limit on the efficiency 
of all thermal engines. Real thermal engines work at 
finite cycle times and lose a finite amount of energy 
due to irreversible cycles and other mechanisms such 
as mechanical friction, heat leak and dissipative pro- 
cesses, etc. Searching for real thermal engines which 
operate with optimal cycles has caught a lot of atten- 
tion. Here "optimal" refers to different optimizations of 
the heat engine, such as maximum efficiency, maximum 
power, maximum entropy production [2] and maximum 
work[3], etc. Of all these optimizations, the efficiency 
of thermal engines at maximum output power is a very 
practical problem and has been extensively studied in 
the literature [UlJ-lal- The efficiency of a quantum ther- 
mal engine operating at maximum power has also re- 
cently been studied 1ZD - 

One of the most important results addressing the ef- 
ficiency of a thermal engine at maximum power was 
given by Curzon and Ahlborn in 19750]. Here we 
briefly review their result first. They made the assump- 
tion that during the time that the working medium is in 
contact with the hot(cold) reservoir, the amount of heat 
exchanged is proportional to the temperature difference 
between the working medium and the reservoirs, and 
also to the time duration of the processes. During the 
heating process, which lasts time t\, the amount of heat 



Wx absorbed by the system is 

W x =Mi(T x -T lw ), 



(1) 



where T\ is the temperature of the heat source, T\ w the 
temperature of the working medium and kx the heat 
transfer coefficient of the heating process. Similarly, for 
the cooling process which lasts time ti, the working 
medium releases heat W2 



W 2 = k 2 t 2 (T 2w -T 2 ). 



(2) 



Here T 2 is the temperature of the cold source, T 2w the 
temperature of the working substance and k 2 the heat 
transfer coefficient of the cooling process. The reversibil- 
ity of the adiabatic stages requires 



Wx W 2 



Tx, 



T 2 



(3) 




This leads to a relationship between tx and t 2 . By maxi- 
mizing the power output of the system, they derived the 
famous Curzon- Ahlborn (CA) formula for the efficiency 
of the thermal engines at maximum power as: 



(4) 



The CA formula describes the thermal engines of power 
plants very well GHH and all the parameters here have 
clear physical meanings. However, as pointed out by 
Ref. |l|], the CA formula is neither exact nor universal, 
and it gives neither an upper bound nor a lower bound. 

In Ref. (10, the authors considered a Carnot thermal 
engine performing finite-time cycles. They assume that 



2 



the amount of heat absorbed by the system per cycle 
from the hot(cold) reservoir is given by 



and 



Q h =T h (AS - — + ...), 

Th 



Q c = T C (-AS - 



(5) 



(6) 



where 7), c is the temperature of the hot(cold) reservoir 
and Th, c the time during which the thermal machine is in 
contact with the hot(cold) reservoir. The second terms of 
Eqs.||5]l and (0 give the extra entropy production per cy- 
cle when the system deviates from the reversible regime. 
By maximizing the power, the efficiency of the engine 
can be derived. The upper and lower bounds of the ther- 
mal efficiency at maximum power are derived when the 
ratio £h/£ c approaches and oo respectively. Espos- 
ito et al's result agrees well with the observed efficien- 
cies of thermal plants [lj, |8y. However, why the work- 
ing medium releases less heat for longer contact times 
with the cold reservoir as indicated by Eq.© was not 
explained. Also, in both Ref.[l] and [4], only Carnot en- 
gines were studied. 

In this paper,we study a more general and realistic 
thermal engine model and derive its efficiency bounds 
at maximum power. In both Ref.[l] and [4], the temper- 
ature of the working medium does not change during 
heat transferring processes which is not true for either a 
realistic system, such as thermal plants, or for other heat 
engine models, such as the Otto, Joule-Brayton,Diesel, 
and Akinson engines[9]. Instead, we simply assume 
that heat transfer by a thermal engine is described by 
Newton's law of cooling, thus it does not have to be 
isothermal anymore. Furthermore, we also take into ac- 
count the fact that the thermal capacities of the working 
medium in realistic systems usually could be quite dif- 
ferent at high and low temperatures! 10]. This is also be 
motivated by heat engines, such as the Diesel and Akin- 
son type, for which the thermal capacities are different 
at the two different thermal stagesl9J]. 

With these two modifications, we argue that our 
model is not only more realistic, but also more general. 
Since the efficiency and its bounds are derived by con- 
sidering heat exchange processes during which the tem- 
perature of the working medium could be close to or far 
away from isothermal, our model could simulate heat 
transferring not only in Carnot engines, as in Ref.[l] 
and [4], but also some other engines such as Otto, Joule- 
Bray ton, Diesel, and Akinson as described in Ref.[9]. 

The organization of the paper is as follows, we will 
first describe Newton's law of cooling and derive the 
corresponding entropy and heat formulas, and then 
study the thermal efficiency at maximum power for two 
limiting cases. 



HEAT TRANSFER AND ENTROPY PRODUCTION 
BASED ON NEWTON'S LAW OF COOLING 

We assume heat transferred by thermal engines in 
contact with a heat source is described by Newton's law 
of cooling: 



(7) 



where c is the heat capacity, rn medium mass, T medium 
temperature, T s heat source temperature, h heat trans- 
fer coefficient, and A contact area. For convenience, we 
denote hA by k. Though Newton's law of cooling is 
quite simple, many other heat transfer laws can be sim- 
plified to it if the temperatures of the objects are high 
while the temperature difference between them is small. 
Based on this assumption, we consider a thermal engine 
working between hot and cold reservoirs at tempera- 
tures Th and T c respectively, and the initial temperature 
of the working medium is Tho(T c o) at the beginning of 
the heating(cooling) stage. The solution to Eq.lO gives 
the temperature of the working medium at time t: 

T(t) = T h + (T h0 -T h )e~^ = T h +(T h0 -T h )e~^ , (8) 

where S/i = Chmjkh- Assuming that the time during 
which the working medium is in contact with the high 
temperature source is 77,, the entropy produced during 
the heating process can be evaluated straightforwardly 
as: 



AS h = - 



h dQ Th - xe 

— — = Chin m 

T(t) T h0 



(9) 



where x = Th — T^o- The heat exchanged between the 
working medium and the high temperature source is 
given by 

Qh= f k h {T h -T h0 )e-^dt = c h mx(l-e-^). (10) 
Jo 

Here and from now on we take the convention that 
Q > means absorbing and Q < releasing heat. Sim- 
ilarly, the entropy production and heat exchange of the 
working medium during the cooling process are given 
by 



AS C = c c m In 



T 



ye 



T C Q 



Q c = -c c my(l-e 



(11) 
(12) 



where y = T c q — T c and E c = c c m/k c . After a thermo- 
dynamic cycle, the system returns to its initial state, and 
the total entropy change of the working medium should 
be zero AS h + AS C = [4], which leads to 



ln[( 



ye 



T c o 



.)C( 



xe 



Tho 



(13) 



3 



By noting that T ha = T h -x, T c0 = T c + y and by defining 
7 = Ch/c c , Eq. ((T3b is reduced to 

V T c + y A r h -a; ; V ' 

The power output and the efficiency of the thermal en- 
gine are given by 



P 



Qh + Qc 

Th + T c 



V = 1 + 7^- 
Hh 



(15) 
(16) 



Generally, the efficiency r\ m at maximum power output 
can be derived using the constraint of Eq.l fT4)l . However, 
Eq. l(l4") is a transcendental equation which can not be 
solved analytically. In what follows we will focus our 
discussions on two special cases. 



EFFICIENCY AND ITS BOUNDS IN TWO SPECIAL 
CASES 

Case I: Long contact time limit: r/E — > oo 
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FIG. 1: (Color online) A comparison of upper and lower 
bounds for the long contact time limit between this work and 
[J] . The black dot-dashed line denotes the CA efficiency. The 
red dashed lines denote the upper and lower bounds of the 
thermal efficiency derived in Ref. ]1]. The orange solid lines 
denote the bounds derived in this paper. 

In this case, the contact time is long enough that the 
working medium can exchange heat sufficiently with 
the reservoirs. Therefore the final temperature of the 
working medium is close to the heat reservoir and 
quite different from its initial temperature. Numerically, 
when r/S - 5, we have \T - T hQ \/\T h - T h0 \ ~ 0.007 



and \T - T c0 \/\T C - T c0 \ ~ 0.007. Thus, when r/S is suf- 
ficiently large, which is supposed to be the case studied 
in Ref.[l], exp (— t/E) can be safely ignored, and Eq. ((T4t 
is reduced to 



(_^£_)(_5i!_)7 = 1. 
K T B + y )K T h -x' 



(17) 



By plugging Eq.l fT7|l into Eq. l(T0) and (12), using 
Eq. l(T5) then the output power is given by 



P = m 



c h x - c c T c [(^y - 1] 



Th + T c 

Let dP/dx = 0, P is maximized when 

x = T h [\- 

J-h 



(18) 



(19) 



Therefore, the efficiency at maximum power is given by 

l r l-£ 



Tjm = 1 



1 



7 l-(§)^ 



--!] = ! — 



Vc 



7 1 - (1 - Vc) 



+7 

(20) 

From the above expression we see that r\ m decreases as 7 
increases. For the symmetric dissipation in which 7 = 1, 
r] m becomes 



Vn 



= 1 - 



(21) 



Interestingly, the CA efficiency is recovered though the 
situation is quite different. Expanding r\ m in series of rj c , 
we have 



Vm = 



Vc 



12 



(1 



1 + 7 



(22) 



The coefficient of the second order term lies between 
1/12 and 1/6, while in Ref.[l] this term is between and 
1/4 which indicates a tighter bound here . The lower 
and upper bounds of rj m in this case are given by 



1 



In (1 - Vc ) 



< Vm < H • (23) 



A comparison of the upper and lower bounds for the 
long contact time limit between this work and results 
derived in Ref.[l] is shown as FIGQ] We see that the 
limits derived here give much tighter bounds than those 
derived in Ref.|l|]. 

We emphasize again that our model does not only ap- 
ply to Carnot engines. Since the final temperature of the 
working medium after heat exchange can be quite dif- 
ferent from its initial temperature, it is not necessarily 
an isothermal process and thus the engine does not need 
to be a Carnot type engine. It can also simulate the en- 
gines described in Ref.|9J], if we take Ch = c c — c v , it is 
the Otto engine, c/, = c c = c p , the Joule-Brayton engine, 



4 



Ch = Cp,c c = c v , the diesel engine and Ch. = c v ,c c = c p , 
the Akinson engine. We can recover all the thermal ef- 
ficiencies at maximum power derived in Ref.|9]. Cor- 
respondingly, the bounds derived in this section should 
apply to those four types of engines mentioned above in 
practical conditions. 



Case II: Short contact time limit: r/E — > 

In this case, the heating and cooling processes are both 
short. Therefore the final temperature of the working 
medium after transferring heat is very close to its initial 
temperature. This is approximately what was studied in 
Ref.[4] where temperature of the working medium does 
not change during heat transfers. Numerically, one can 
estimate that if t/E - 0.1, then \T - T ha \/\T h - T h0 \ ~ 
0.9 and \T - T cQ \/\T c - T c0 \ ~ 0.9. We solve Eq.Qll by 
expanding it as series of the infinitesimal variable (r/E) 
and matching both sides of the equation order by order 
(We always keep the same order of t/j/E/i and r c /E c ). 
As t/E — > 0, to the zeroth order of r/E, Eq. (TT4]) simply 
gives 1 = 1 which is trivial. To the first order of t/E, 
Eq-O gives 



72; 



To obatian the expression for the power output P, plug 
Eqs.ll2"7t and 11281 into Eq. lfl5] |, and expand the expres- 
sions of Qh and Q c to the first order of t/E again, we 
have 



P = 



£ j_ (T h +~fx)Y; c 



(29) 



T h -x 



P is maximized by letting dP/dx = 0. Note < x < Th 
and x < ^Tc, the unique allowed solution of x is given 
by 



(T e S fe +T fe S e )(l+ 7 ) _ ^j,^ 



(T h +7T c )(£ h +E c ) 



S h +S c 



(30) 



Therefore the thermal efficiency at maximum power is 
given by 



7Sc-S„ 1 / (r e S h +T h S e )(l+7) 
T c S„+S c 1 + (T h +7T C )(S, 1+ S c ) 
Vm = 1 — — ; ' (31) 

T h 7 Sc-s fe _j_ 1 _ / (r e s h +r h s c )(i+7) 

S h +S c "f" 1 V (T, 1 + 7 T c )(S h + S <? ) 



(24) By defining (3 = kh/k ci and using the relation r\ c = 1 



Now the amounts of heat exchanged by the system dur- 
ing the heating and cooling processes are given by 



Q h = c h mx—, 
Q c = -c c my— 



(25) 
(26) 



The above equations agree with the fundamental equa- 
tions listed at the beginning of Ref. J3]. Thus if we 
continue our straightforward calculation, we simply re- 
cover the same results in Ref. [4], including the CA effi- 
ciency. Now, we continue to expand Eq. l(T4)l to the sec- 
ond order of t/E, we obtain another simple relation 



V 
x 



(27) 



Combining this with Eq. (|24l l, we get 

Tc_ _ T h + JX T h 

E c T h - x T,h 




0.2 0.4 0.6 0.8 1 
\ 

FIG. 2: (Color online) A comparison of Eq.ll35l (red solid line) 
with CA efficiency (black dot-dashed line). 



(28) Tc/T h , v m can be expressed as 



Vm = 1 - 7(1 - Vc) 



+ 7(1 - vcW + 7) - VI 1 + 7(1 - VcW + 7) 



+ 7(1 " Vc)](l + 7) - \/W + 7(1 " VcW + 7) ' 



(32) 



5 



= 1 - 



2 72 



21 
7l' 



(1 - n' c ) - 1 + ^^(l + 7i 2 (l-^) 2 )(7i 2 +7 2 4 ) 



7l + 7 1 2 (l-^) 



(33) 



Moreover, r\ m can be expanded in a series of r\ c as 



1 1 1 _J- 

2 Vc+ 8 ( l + 7 //3 + 1 + 7 



+ T-rrX + 0(vi). (34) 



The coefficient of the first order term of r\ m is 1/2, and 
the coefficient of the second order term lies in the range 
between and 1/4. In the symmetric case where /3 = 1 
and 7=1, we have 



7?m 



?7c(2 - r] c ) 
4 - 3t7 c 



(35) 



When expanding as a series in rj c , the coefficient of the 
second order term is 1/8. Those results agree with the 
expansion of CA efficiency jj]. To the third order of r/ c , 
the difference of r] m from tjca is ?/c/32 + 0{rQ). A com- 
parison between rj m and the CA efficiency and our result 
for the symmetric case is shown in FIGf2J 

Now, we estimate the bounds of r\ m . In the limits 7 = 
or 7 = 00 while (3 is finite, we recover the lower and 
upper limits of r) m given in Ref. (J. 



Vc , 

"2 " % 



< 



(36) 



Interestingly our bounds on are obtained in the short 
contact time limit (r — > 0) while the same results were 
obtained in the long contact time limit (t — > 00) in Ref. 

It is easy to verify that rj m decreases as 7 increases, 
but increases as j3 increases. This means the larger the 
ratio between heat capacities of the working medium at 
the hot and cold reservoirs, the lower the efficiency at 
maximum output power. 



CONCLUSIONS 

In summary, we presented an analysis of thermal ef- 
ficiency and its bounds at maximum power for thermal 
engines for which the heat transferring processes are de- 
scribed by Newton's law of cooling. In the long con- 
tact time limit, CA efficiency is recovered for symmetric 
thermal capacity and two tighter bounds on the ther- 
mal efficiency are derived. The model can simulate Otto, 



Joule Brayton, Diesel and Atkinson engines in the long 
contact time limit. In the short contact time limit, we 
recover the famous CA efficiency in the first order cal- 
culation. When we proceed to the second order calcu- 
lation, we derived a different efficiency formula and re- 
covered the efficiency bounds at maximum power given 
by Espositi, et al. In both limits, the thermal efficiency 
is found to decrease as 7 = Ch/c c increases. This might 
be helpful for choosing a suitable working medium and 
working temperatures when designing a thermal en- 
gine whose heat transfer can be approximated by New- 
ton's law of cooling. Other cases such as those associ- 
ated with intermediate thermal contact time and differ- 
ent heat transfer laws are being investigated further. 
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